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TITLE 

CONTROL METHOD AND SYSTEM FOR MOTOR 
BACKGROUND OF THE INVENTION 

1. Field of the invention 

The present invention relates to a control method and system for a 
motor, more particularly, to a control method and system using a micro- 
programming & memory controller for a motor. 

2. Description of the Related Art 

Fig. 1 shows a conventional PI-PWM (Proportion Integration-Pulse 
Width Modulation) current control system 11 for a motor 10. The 
conventional PI-PWM current control system 11 comprises: a first 
difference operator 12, a PI controller 14, a second difference operator 15, 
a PWM comparator 16 and an inverter 18. The control system 11 receives 
three phase reference currents 20 to the first difference operator 12, and 
detects three phase actual currents 22 from a motor 10. The first difference 
operator 12 calculates three phase difference currents between the three 
phase reference currents 20 and three phase actual currents 22. 

The three phase difference currents are amplified by the PI controller 

14 to obtain three phase control voltages. The second difference operator 

15 calculates difference voltages between the three phase control voltages 
and three phase triangle-wave signals 24. The PWM comparator 16 
receives the three phase difference voltages and outputs three phase switch 
signals to the inverter 18. According to the three phase switch signals, the 
inverter 18 outputs three phase input voltages to the motor 10 and controls 
the motor 10. 

However, the conventional PI-PWM current control system 1 1 has 
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the current phase-lag problem. Referring to Figs. 2A, 2B, 3A, 3B and 3C, 
they show the comparison between the three phase actual currents and the 
three phase reference currents, according to the conventional PI-PWM 
current control system. Fig. 2A shows A phase reference current I al> refer . 
Fig. 2B shows A phase actual current I al>actual . Fig. 3 A shows a comparison 
between A phase reference current I al>refer and A phase actual current 
V actual • Fig. 3B shows a comparison between B phase reference current 
Ibi.refer and B phase actual current I b i, ac tuai • Fig. 3C shows a comparison 
between C phase reference current I cl)refer and C phase actual current I cUctuaI 
. As shown in Figs. 2A, 2B, 3A, 3B and 3C, the current phase-lag problem 
exists between the three phase reference currents I al , refer , I bl , refer , I cl refer and 
three phase actual currents I al , actual , I b];actual , I cl actual . 

Besides, when a software is utilized to perform the operation of the 
PI controller 14, the conventional PI-PWM current control system 1 1 has 
the high noise problem. Figs. 4A, 4B, 5A, 5B and 5C show the comparison 
between the three phase actual currents and the three phase reference 
currents, according to the conventional PI-PWM current control system 1 1 
using a software to perform the PI controller 14. Fig. 4A shows A phase 
reference current I a2 , refer , Fig. 4B shows A phase actual current I a2 , actual , 
Fig. 5A shows a comparison between A phase reference current I a2!refer and 
A phase actual current I a2 , actua , , Fig. 5B shows a comparison between B 
phase reference current I b2 refer and B phase actual current I b2 , actual , Fig. 5C 
shows a comparison between C phase reference current I c2>refer and C phase 
actual current I c2 , actua i . As shown in Figs. 4A, 4B, 5A, 5B and 5C, the three 
phase actual currents I al)actual , I bl actual , I cl actual have the high noise problem. 

Therefore, it is necessary to provide an innovative and progressive 
control method and system to solve the above problem. 

SUMMARY OF THE INVENTION 

One objective of the present invention is to provide a control method 
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and system for a motor. The control method comprises the steps of: (a) 
determining a motor parameter; (b) providing a d-axis reference current and 
a q-axis reference current; (c) detecting actual currents of the motor and 
converting to a d-axis actual current and a q-axis actual current; (d) 
calculating a d-axis counter electromotive force and a q-axis counter 
electromotive force at a sample period according to the motor parameter, 
the d-axis actual current and the q-axis actual current at the sample period, 
the d-axis actual current and the q-axis actual current at a last sample 
period, and a d-axis voltage and a q-axis voltage at the last sample period; 
and (e) calculating a d-axis voltage and a q-axis voltage at the sample 
period according to the motor parameter, the d-axis actual current and the 
q-axis actual current at the sample period, the d-axis reference current and 
the q-axis reference current at a next sample period, the d-axis counter 
electromotive force and the q-axis counter electromotive force calculated at 
the sample period. 

According to the control method of the invention, the d-axis voltage 
and the q-axis voltage are calculated and converted to three phase control 
voltages. The three phase control voltages are processed to control the 
motor. The control method of the invention not only preserves the merit of 
constant switching frequency of the conventional PI-PWM control system, 
but also can eliminate the motor current phase-lag problem of the 
conventional PI-PWM control system. Also the control method of the 
invention will have the following characteristics and economic advantages 
such as high accuracy, fast response, low cost and robustness etc. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 shows a block diagram of a conventional PI-PWM current 
control system for a motor. 

Fig. 2 A shows A phase reference current I a i, refer waveform. 

Fig. 2B shows A phase actual current I a i, actual waveform, according to 
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the conventional PI-PWM current control system. 

Figs. 3A to 3C show comparison waveforms between the three phase 
reference currents and the three phase actual current, according to the 
conventional PI-PWM current control system. 

Fig. 4A shows A phase reference current I a2 , refer waveform. 

Fig. 4B shows A phase actual current I a 2, ac tuai waveform, according to 
the conventional PI-PWM current control system using a software to 
perform the operation of PI controller. 

Figs. 5A to 5C show comparison waveforms between the three phase 
reference currents and the three phase actual current, according to the 
conventional PI-PWM current control system using a software to perform 
the operation of PI controller. 

Fig. 6A shows a block diagram of a control system for a motor, 
according to the invention. 

Fig. 6B shows a block diagram of a micro-programming & memory 
controller, according to the invention. 

Fig. 7 shows a flow chart of a control method for a motor, according 
to the invention. 

Fig. 8 A shows A phase reference current I a3> re f er waveform. 

Fig. 8B shows A phase actual current I a 3, act uai waveform, according to 
the control system of the invention using a correct leakage inductance. 

Figs. 9 A to 9C show comparison waveforms between the three phase 
reference currents and the three phase actual current, according to the 
control system of the invention using a correct leakage inductance. 

Fig. 10A shows A phase reference current I a4j refer waveform. 
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Fig. 1 OB shows A phase actual current I a 4,actuai waveform, according 
to the control system of the invention using an excessive leakage 
inductance. 

Figs. 11A to 11C show comparison waveforms between the three 
phase reference currents and the three phase actual current, according to the 
control system of the invention using an excessive leakage inductance. 

Fig. 12A shows A phase reference current I a5) refer waveform. 

Fig. 12B shows A phase actual current I a 5, ac tuai waveform, according 
to the control system of the invention using an inadequate leakage 
inductance. 

Figs. 13 A to 13C show comparison waveforms between the three 
phase reference currents and the three phase actual current, according to the 
control system of the invention using an inadequate leakage inductance. 

Fig. 14A to 14C show three phase actual currents, according to the 
control system of the invention using a software to tune the leakage 
inductance. 

Fig. 14D shows the change of the leakage inductance, according to 
the control system of the invention using a software to tune the leakage 
inductance. 

DETAILED DESCRIPTION OF THE INVENTION 

The control method of the invention can be utilized to control an AC 
motor, more particularly to control an induction motor so that three phase 
actual currents are equal to three phase reference currents. In the 
embodiment, the induction motor is taken as an example to describe the 
operation of the control method and system of the invention. However, the 
control method and system of the invention are not limited to control the 
induction motor, the other AC motor and DC motor can be controlled by 
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the control method and system of the invention. 

Firstly, the induction motor dynamic equations on the synchronous 
rotating d-q frame are shown as follows Equation (1). 
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wherein L a is defined as a leakage induction. 
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L ss = ^l_ji M.=leakage inductance 



Equation (1) is re-arranged, and a d-axis voltage V ds and a q-axis 
voltage V qs are obtained as follows. 
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(3) 



A d-axis counter electromotive force E ds and a q-axis counter 
electromotive force E qs are defined as follows. 



E ds =-<o r ^\ r + 



E qs =-HO r ^ dr+ -i qs -r r ^V 



(4) 
(5) 
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Therefore, Equation (2) and (3) can be simplified as follows. 



Vqs=L °lf + LoC ° cids+E * ( 7 ) 



Wherein di ds /dt and di qs /dt can be calculated by a d-axis reference 
current i dS(I . ef , a q-axis reference current i qs ref , a d-axis actual current i dS) actuai 
and a q-axis actual current i qS; actua i as follows. 



di* _ i ds,ref( t + T s)- i ds,actual( t ) _ i^error ( 8 ) 

dt T T 

s s 

digs _ l qs,ref(t + ^s) actual^) = ^gs2 t error ( 9 ) 



dt T 



s 



Wherein T s is a sample time, and T s =0.2ms. The i dS;actua i (t) means a 
ids,actuai value at a sample period. The i qS;actua i (t) means a i qs , ac tuai value at a 
sample period. The i ds> ref (t+T s ) means a i^f value at a next sample 
period. The i qS)ref (t+T s ) means a i qs ref value at a next sample period. 
Therefore, di ds /dt and di qs /dt can be calculated easily by Equations (8) and 
(9). 

The d-axis counter electromotive force E ds and the q-axis counter 
electromotive force E qs can be derived from Equations (6) and (7) as 
follows. 



Eds =V ds -L 0 -^ + L CT a> e i qs (10) 
di„ 



E qs = V qs-L a -^ L -L (J CD e i ds 



(11) 
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Wherein di ds /dt and di qs /dt can be calculated by the d-axis actual current 
ids, actual and the q-axis actual current i qSj actual as follows. 

di ds _ ids, actual (0 ~ hs,actual (? ~ ) _ ^dsLerror 



(12) 

^qs _ ^qs, actual (0 ~ ^qs, actual ~ ) _ ^qs\, error (13) 



a t; r 



Wherein T s is a sample time, and T s =0.2ms. The id S ,actuai (t) means a 
ids,actuat value at a sample period. The i qs , a ctuai (X) means a i qs ,actuai value at a 
sample period. The i ds , actU ai (t-T s ) means a id S ,actuai value at a last sample 
period. The i qs , a ctuai (t-T s ) means a i qs , a ctu a i value at a last sample period. 
Therefore, di^/dt and di qs /dt in Equations (10) and (11) can be calculated 
easily by Equations (12) and (13). 

Therefore, the d-axis voltage V ds and the q-axis voltage V qs can be 
calculated by Equations (6) to (13) according to the leakage inductance, the 
d-axis reference current i^f , the q-axis reference current i qSjre f , the d-axis 
actual current id S ,actuai and the q-axis actual current i qS)act uai • 

Referring to Fig. 6A, Fig. 6B and Fig. 7, the control system 60 
comprises: a micro-programming & memory controller 61, a PWM 
comparator 62 and a switch-mode inverter 63. The micro-programming & 
memory controller 61 comprises: an input means 611, a detecting means 
612, a first calculating means 613, a second calculating means 614, a 
converting means 615 and a third calculating means 616. The input means 
61 1 is used for receiving a motor parameter, a d-axis reference current i ds , r ef 
and a q-axis reference current i qs , re f- The motor parameter can be a leakage 

( L 2 ^ L L -L 2 

inductance L a , and the expression is l o = l - _ M |l cc = ss n The 



L L L 

leakage inductance L a can be a constant. As shown in step 701 of Fig. 7, 
the motor parameter is determined. 

The input means 611 can receive three phase reference currents 601, 
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then the three phase reference current are converted to the d-axis reference 
current and the q-axis reference current i qs>re f , as shown in step 702 of 
Fig. 7. Then, it is easy to obtained the d-axis reference current i^f (t+T s ) 
and the q-axis reference current i qSjre f (t+T s ) at a next period in Equations (8) 
and (9). 

The detecting means 612 is used for detecting three phase actual 
currents 602 from a motor 64 and converting the three phase actual currents 
to a d-axis actual current i ds , ac tuai and a q-axis actual current i qs , a ctuai , as 
shown in step 703 of Fig. 7. At a sample time, the d-axis actual current 
ids,actuai (t) and the q-axis actual current 

Aqs^ctual 

(t) can be obtained in real- 
time, and at a last sample time (t-T s ) the values of the d-axis actual current 
ids,actuai (t-T s ) and the q-axis actual current i qs , a ctuai (t-T s ) are stored in a 
memory. 

Therefore, the di ds /dt and di qs /dt in Equations (12) and (13) can be 
calculated. The id S ,actuai (t) is the id S ,actuai value at a sample period, the 
i qs ,actuai(t) is the i qs , act uai value at a sample period, and the i dSj actU ai (t-T s ) is the 
stored ids,actuai value at a last sample period. The i qs , ac tuai (t-T s ) is the stored 
i qs ,actuai value at a last sample period. 

The first calculating means 613 is utilized to calculate a d-axis counter 
electromotive force E ds and a q-axis counter electromotive force E qs at a 
sample period according to the motor parameter L a , the d-axis actual 

current i ds , a ctuai (0 and the q-axis actual current i q s,actuai 

(t) at the sample 

period, the d-axis actual current i ds , ac tuai (t-T s ) and the q-axis actual current 
iqs, actual (t-T s ) at a last sample period, and a d-axis voltage V ds and a q-axis 
voltage V qs at the last sample period. Therefore, the d-axis counter 
electromotive force E ds and the q-axis counter electromotive force E qs are 
calculated by Equations (10) and (11). The method for calculating the d- 
axis counter electromotive force E ds and the q-axis counter electromotive 
force E qs is shown in step 704 of Fig. 7. 
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The second calculating means 614 is used for calculating a d-axis 
voltage Vds and a q-axis voltage V qs at the sample period according to the 
motor parameter L a , the d-axis actual current id S ,actuai (t) and the q-axis 
actual current i qs , a ctuai (t) at the sample period, the d-axis reference current 
ids,ref (t+T s ) and the q-axis reference current i qs , r ef ( t+ T s ) at a next sample 
period, the d-axis counter electromotive force and the q-axis counter 
electromotive force E qs at the sample period calculated by the first 
calculating means 612. Therefore, the d-axis voltage and the q-axis 
voltage V qs are calculated by Equations (6) and (7). The method for 
calculating the d-axis voltage V ds and the q-axis voltage V qs is shown in 
step 705 of Fig. 7. The d-axis voltage V ds and the q-axis voltage V qs can be 
used to control the motor 64. 

At the first sample period, the d-axis counter electromotive force E ds 
and the q-axis counter electromotive force E qs are equal to zero. Then, the 
d-axis voltage V ds and the q-axis voltage V qs can be calculated by 
Equations (6) and (7) at the first sample period using E ds =0 and E qs =0. At 
the second sample period, the d-axis counter electromotive force E ds and 
the q-axis counter electromotive force E qs can be calculated by Equations 
(10) and (11) using the d-axis voltage V ds and the q-axis voltage V qs at the 
first sample period. Then, the d-axis voltage V ds and the q-axis voltage V qs 
can be calculated by Equations (6) and (7) at the second sample period 
using the d-axis counter electromotive force E ds and the q-axis counter 
electromotive force E qs calculated at the second sample period. The above 
steps are repeated to calculate the d-axis voltage V ds and the q-axis voltage 
V qs every sample period so as to control the motor 64 of Fig. 6A. 

The converting means 615 is used for converting the d-axis voltage V ds 
and the q-axis voltage V qs to three phase voltage V as ,V bs , V cs , as shown in 
step 706 of Fig. 7. The three phase voltage are converted by the third 
calculating means 616 to obtain three phase control voltage V a)Contro i , 
V b ,controi and V c>contro i > as shown in step 707 of Fig. 7, according to the three 
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phase voltages V as ,V bs , V cs , a peak voltage of a comparing voltage V Tri _p 
and a DC voltage of the switch-mode inverter V d as follows. 



The three phase control voltages are compared with a comparing 
voltage 603 by the PWM comparator 62, as shown in step 708 of Fig. 7. 
The comparing voltage 603 is a triangle wave voltage signal. After the 
PWM comparator 62, three phase switching signals are obtained to input 
and control the switch-mode inverter 63, as shown in step 709 of Fig. 7. 
The switch-mode inverter 63 outputs three phase input voltages to the 
motor 64, as shown in step 710 of Fig. 7. 

Therefore, the three phase input voltages of motor can be calculated by 
the control method and the control system, according to the invention. The 
motor 64 can be controlled by the three phase input voltages so as to 
perform the objective of the three phase actual currents being equal to the 
three phase reference current. 

Figs. 8A, 8B, 9A, 9B and 9C show the comparison between the three 
phase actual currents and the three phase reference currents, according to 
the control system 60 and the control method of the invention. Fig. 8A 
shows A phase reference current I a3> re f er > Fig. 8B shows A phase actual 
current I a 3,actuai > Fig- 9 A shows a comparison between A phase reference 
current I a 3, refer and A phase actual current I a 3, actual > Fig. 9B shows a 
comparison between B phase reference current prefer and B phase actual 
current factual , and Fig. 9C shows a comparison between C phase 
reference current I C 3, re fer and C phase actual current I C 3, a ctuai • Figs. 8A, 8B, 
9A, 9B and 9C are obtained by utilizing the Matlab software to perform the 
operation of the control system and using a correct leakage inductance. As 
shown in Figs. 9A, 9B and 9C, the three phase actual currents I a3 ,actuai , 
factual , I c 3,actuai are almost equal to the three phase reference currents I a 3,refer 




1 4 
1 5 
1 6 
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. Ib3,refer , I C 3,refer, respectively. Therefore, the control system can solve the 
current phase-lag problem in the conventional PI-PWM current control 
system. 

To prove the robustness of the control system according to the 
invention, the Matlab software is also utilized to perform the operation of 
the control system using an excessive value of leakage induction and an 
inadequate value of leakage induction. 

Figs. 10A, 10B, 11 A, 11B and 11C show the comparison between the 
three phase actual currents and the three phase reference currents, 
according to the control system 60 and the control method of the invention 
and using an excessive value of leakage induction 1.5 L a . Fig. 10A shows 
A phase reference current I a4 , re f er , Fig. 10B shows A phase actual current 
factual , Fig. 1 1 A shows a comparison between A phase reference current 
Vrefer and A phase actual current I a4)actual , Fig. 11B shows a comparison 
between B phase reference current I M>refer and B phase actual current I b4 , actuaI 
, and Fig. 11C shows a comparison between C phase reference current 
Urefer and C phase actual current I c4)actual . As shown in Figs. 1 1 A, 1 IB and 
11C, even using the excessive leakage induction 1.5 L ff , three phase actual 
currents I a4jactU ai , factual , I C 4,actuai are still almost equal to the three phase 
reference currents I a4 , refer , I b4jrefer , l c4 refeo respectively. 

Figs. 12A, 12B, 13A, 13B and 13C show the comparison between the 
three phase actual currents and the three phase reference currents, 
according to the control system 60 and the control method of the invention 
and using an inadequate value of leakage induction 0.5 L a . Fig. 12A 
shows A phase reference current I a5 , refer , Fig. 12B shows A phase actual 
current I a5)actua , , Fig. 13 A shows a comparison between A phase reference 
current I a5>refer and A phase actual current I a5>actuaI , Fig. 13B shows a 
comparison between B phase reference current I b5>refer and B phase actual 
current I b5) actuai , and Fig. 13C shows a comparison between C phase 
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reference current I c5>refer and C phase actual current I c5;actual . As shown in 
Figs. 13 A, 13B and 13C, even using the inadequate leakage induction 0.5 
L a , three phase actual currents I a5 , actual , I b5actuaI , i c5 actua) are still almost 

equal to the three phase reference currents I a5 , refer , I b5jrefer , Wferj 
respectively. 

Therefore, as shown in Fig. 8 A to Fig. 13C, within ±50% error of the 
leakage inductance L a , the three phase actual currents can be almost equal 
to the three phase reference currents. The control method of the invention 
not only preserves the merit of constant switching frequency of the 
conventional PI-PWM control system, but also can eliminate the motor 
current phase-lag problem of the conventional PI-PWM control system. 
Also, the control method of the invention will have the following 
characteristics and economic advantages such as high accuracy, fast 
response, low cost and robustness etc. 

In the above embodiments, the leakage inductance L a is the only 
variable. The leakage inductance L a is determined usually by the 
parameter test of the motor. If the parameter test is not correct, the value of 
the leakage inductance L a may be above the correct value in the error 
50%. If the incorrect leakage inductance L a with above ±50% error is 
utilized to calculate the d-axis voltage and the q-axis voltage V qs , the 
three phase actual currents may be not equal to the three phase reference 
currents. Therefore, the leakage inductance L a must be real-time auto- 
tuning as follows. 



At a sample period, Equations (6) and (7) are modified as follows. 

Vdsl = L ° if" " Lo<0 ' [ ^ + E <*' = m <» L ° + E dsi ( j 7 ) 
=L -- d f + L « a) ei dsI +E qsl =m ql L 0 +E qsl (18) 
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Wherein in Equations (17) and (18) the variables are the leakage 
inductance L a , E ds] and E qs , , the and n^i are defined and obtained as 

follows = (^i- -CO i , ) . fn = (^fL + aJ i \ 

dl v e>l ' dt e'dsi > • 

At a next sample period, Equations (6) and (7) are modified as 
follows. 

m 42 L o +E*2 (19) 

m q2 L a +E qs2 (2 0) 



V*2=L <r -^-L a a> e i qsl+ E tb2 = 



Vqs2 — L a 



dt 
dt 



qsl 



Wherein in Equations (19) and (20) the variables are the leakage 
inductance L a , E ds2 and E qs2 , the m^ and m^ are defined and obtained as 

follows m = (-^2. _ & i \ m t^ivl. + (oi \ 

dt e «' 2 ' ?2 v j t ^ UJ e l ds2>- 

Equation (19) minus Equation (17) is as follows. 

( v ds2 ~V dsl ) = (m d2 - m dl )L„ + (E M - E*. ) s (m d2 - m dl )L C (21) 

In Equation (21), because the sample time is very small, the 
difference between the E ds2 and E^ can be ignored so that s E w . 
Therefore, the leakage inductance L a d is estimated using d-axis Equations 
(17) and (19). The expression of the leakage inductance L a d is as follows. 

L cd = ~ V*, )/(m d2 - m dI ) (22) 

Similarly, Equation (20) minus Equation (18), and the difference 
between the E qs2 and E qs i can be ignored. The leakage inductance L a is 

estimated using q-axis Equations (18) and (20). The expression of the 
leakage inductance L a q is as follows. 
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Lc, q =(V qs2 -V qsl )/(m q2 ~m ql ) (2 3) 

Therefore, the leakage inductance L a is defined as an average of the 
leakage inductance L ad and the leakage inductance L a q as follows. 

L CT =0.5*(L a , d +L CT J (2 4) 

The leakage inductance L a is calculated by the above Equations 
(17) to (24). The calculated leakage inductance L a is inputted to 
Equations (6), (7), (10) and (1 1) to calculate V ds , V qs , E ds and E qs for two 
sample periods as V ds i , V qs i , V<js2 and V qs2 . Therefore, the value of , 
V qsl , V ds2 and V qs2 in Equations (17), (18), (19) and (20) can be obtained 
to calculated a new leakage inductance L a . Every two sample periods , the 

leakage inductance L a is calculated and updated to input to Equations (6), 

(7), (10) and (1 1) to calculate V ds , V qs , E ds and E qs . Therefore, Equations 
(17) to (24) is calculated every two sample periods to obtained a new 
leakage inductance L a so that the leakage inductance L a can approach the 

correct value. 

Referring to Fig. 14A to Fig. 14D, it shows a simulation using the 
Matlab software. As shown in Fig. 14D, the leakage inductance L a can 

approach to the correct value and be stable within 0.25 second. 

Referring to Fig. 6B, the micro-programming & memory controller 
61 further comprises an fourth calculating means 617 for estimating the 
leakage inductance L a according to the above method. Therefore, by 

using the estimation of the leakage inductance L a the control system 60 

can improve the robustness. 

While an embodiment of the present invention has been illustrated 
and described, various modifications and improvements can be made by 
those skilled in the art. The embodiment of the present invention is 
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therefore described in an illustrative but not restrictive sense. It is intended 
that the present invention may not be limited to the particular forms as 
illustrated, and that all modifications which maintain the spirit and scope of 
the present invention are within the scope as defined in the appended 
claims. 
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